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Cic ABSTRACT 
I 

2 Aims. An unbiased and detailed characterization of the galaxy luminosity function (LF) is a basic requirement in many astrophysical issues: it is 

of particular interest in assessing the role of the environment in the evolution of the LF of different galaxy types. 
^ Methods. We studied the evolution in the B band LF to redshift z ~ 1 in the zCOSMOS 10k sample, for which both accurate galaxy classifications 
I ^ I (spectrophotometric and morphological) and a detailed description of the local density field are available. 

Results. The global B band LF exhibits a brightening of ~ 0.7 mag in M' from z ~ 0.2 toz ~ 0.9. At low redshifts (z < 0.35), spectrophotometric 
late types dominate at faint magnitudes (M^^g > -20), while the bright end is populated mainly by spectrophotometric early types. At higher 
^ redshift, spectrophotometric late-type galaxies evolve significantly and, at redshift z ~ l,the contribution from the various types to the bright end 
of the LF is comparable. The evolution for spectrophotometric early-type galaxies is in both luminosity and normalization: M' brightens by ~ 0.6 
J — mag but (f)' decreases by a factor ~ 1.7 between the first and the last redshift bin. A similar behaviour is exhibited by spectrophotometric late-type 
galaxies, but with an opposite trend for the normalization: a brightening of ~ 0.5 mag is present in M*, while 0* increases by a factor ~ 1.8. 
Studying the role of the environment, we find that the global LF of galaxies in overdense regions has always a brighter M* and a flatter slope. 
• In low density environments, the main contribution to the LF is from blue galaxies, while for high density environments there is an important 
contribution from red galaxies to the bright end. 

The differences between the global LF in the two environments are not due to only a difference in the relative numbers of red and blue galaxies, but 

also to their relative luminosity distributions: the value of M* for both types in underdense regions is always fainter than in overdense environments. 

These results indicate that galaxies of the same type in different environments have different properties. 
^ We also detect a differential evolution in blue galaxies in different environments: the evolution in their LF is similar in underdense and overdense 
. i-H regions between z ~ 0.25 and z ~ 0.55, and is mainly in luminosity. In contrast, between z ~ 0.55 and z ~ 0.85 there is little luminosity evolution 

but there is significant evolution in 0*, that is, however, different between the two environments: in overdense regions 0* increases by a factor 
5^ ~ 1.6, while in underdense regions this increase reaches a factor ~ 2.8. Analyzing the blue galaxy population in more detail, we find that this 

evolution is driven mainly by the bluest types. 

Conclusions. The "specular" evolution of late- and early-type galaxies is consistent with a scenario where a part of blue galaxies is transformed in 
red galaxies with increasing cosmic time, without significant changes in the fraction of intermediate-type galaxies. The bulk of this tranformation 
in overdense regions probably happened before z ~ 1, while it is still ongoing at lower redshifts in underdense environments. 

Key words. Galaxies: evolution - Galaxies: luminosity function, mass function - Galaxies: statistics - Surveys 



1. Introduction 

The COSMOS project (Scoville et al. l2007] l aims to identify the 
physical processes that drive the evolution of galaxies. For exam- 
ple, dynamical processes are likely to play a major role in deter- 
mining the galaxy morphology, whereas dissipative phenomena 
affect the gas content and therefore the star formation, altering 
the galaxy spectral energy distribution. The interplay between 
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these two types of processes is not yet completely understood, 
as well as the role of the environment in the evolution of galaxy 
properties. 

From the observational point of view, the investigation of 
these topics requires: a) high quality images, for deriving ac- 
curate morphological classifications; b) complete multiwave- 
length coverage, for determining spectral energy distributions; c) 
galaxy spectra, for obtaining spectroscopic redshifts (necessary 
for a precise environment description) and measuring spectral 
features (to be used as diagnostics of the gas and stellar physics). 
The COSMOS multiwavelength imaging project provides data 



2 



E.Zucca et al.: The zCOSMOS luminosity function 



for points a) and b), while the zCOSMOS survey (Lilly et al. 
I2007I I was developed to fullfill point c). The combination of 
these data allows us to study with unprecedented completeness 
the properties and evolutionary histories of galaxies as a function 
of their type and environment at the same time. The luminosity 
function is the first and most direct estimator in quantifying this 
evolution. 

Several works have already demonstrated that the global 
luminosity function evolves. The Canadian Network for 
Observational Cosmology field galaxy redshift survey (CNOC- 
2, Lin et al. [19991 ) and the ESO Sculptor Survey (ESS, deLappar- 
ent et al.'2003') derived the luminosity function up to z ~ 0.5 us- 
ing ~ 2000 and ~ 600 redshifts, respectively. At higher redshift, 
the Canada France Redshift Survey (CFRS, LiUy et al. [T995] l 
allowed to study the luminosity function up to z ~ 1.1 with a 
sample of ~ 600 redshifts. 

A major improvement in this field was obtained with the VIMOS 
VLT Deep Survey (VVDS, Le Fevre et al. 2005 ), which detected 
a significant brightening of the M* parameter, amounting to ~ 2 
mag in the B band from z ~ to z ~ 2 (Ilbert et al. l20051 l. using a 
sample of ~ 11 000 spectra. With the same sample, Zucca et al. 
( |2006t studied the evolution in the luminosity function for dif- 
ferent spectrophotometric types, finding a strong type-dependent 
evolution and identifying the latest types as being responsible for 
most of the evolution in the global luminosity function. 
Establishing the role of the environment is more difficult, be- 
cause accurate redshifts and surveys with high sampling rates 
are necessary for reliable density estimates. These constraints 
are satisfied by large local surveys, such as the two-degree Field 
Galaxy Redshift Survey (2dFGRS, Colless et al. 2001 ) and the 
Sloan Digital Sky Surv ey (SD SS, York et al. 2000). From the 
2dFGRS, Croton et al. ( I2005I I measured the dependence of the 
luminosity function on the density contrast (defined in spheres 
of radius 8 h"' Mpc) and galaxy types, finding that the void re- 
gions are dominated by faint late-type galaxies and that the clus- 
ter regions exhibit an excess of very bright early-type galaxies. 
The parameter M* brightens in overdense regions for all galaxy 
types, while the slope a steepens with increasing density for late- 
type galaxies and remains constant for early-type galaxies. 
At high redshift, the situation becomes more complicated: many 
surveys are based on photometric redshifts, which do not allow 
an accurate reconstruction of the density field, and/or lack galaxy 
classifications (morphological and/or spectrophotometric). 
Deep surveys based on spectroscopic redshifts, such as the 
VIMOS-VLT Deep Survey (VVDS, Le Fevre e t al. |2005T l and 
the DEEP2 Galaxy Redshift Survey (Davis et al.|2003l), are able 
to study galaxy properties as a function of the environment to 
z ~ 1.5. Cucciati et al. (2006), using VVDS data, and Cooper 
et al. (I2P06 2007), using DEEP2 data, claimed that the color- 
density relation significantly evolves with redshift, and Cucciati 
et al. (12006 ) found that this relation is established at higher red- 
shift for brighter galaxies. To understand these results, it is nec- 
essary to study the evolution in the luminosity function of differ- 
ent galaxy types in different environments. Until now, the data 
required to complete this study as a function of galaxy type, en- 
vironment, and redshift, has been unavailable. 
In this paper, we study the effect of the environment on the evo- 
lution in the luminosity function of different galaxy types for the 
zCOSMOS 10k bright sample. Parallel studies discuss the evo- 
lution of the luminosity density (Tresse et al. 120091 ), the mass 
function for different types (Pozzetti et al. 12009 ) and in various 
environments (Bolzonella et al. 2009|, the spectrophotometric 
properties as a function of the environment (Cucciati et al. i2009[ 
lovino et al. ,2009) and the morphology-density relation (Tasca 



et al. 120091 Kovac et al. |2009bl ). 

The paper is organized as follows: in Sect. |2| we present the 
data on which this work is based and in Sect. |3]we describe the 
method used to estimate luminosity functions. The results con- 
cerning the luminosity function evolution are described in Sect. 
|4|and the contribution of the different galaxy types is presented 
in Sect. |5] paying particular attention to early-type galaxies in 
Sect. |6] The role of the environment is discussed in Sect. |7| and 
the results are summarized in Sect. |8] 

Throughout this paper we adopt a flat Q,„ = 0.25 and Qa = 0.75 
cosmology, with Ho = 70 km s"' Mpc"'. Magnitudes are given 
in the AB system. 



2. Data 

The zCOSMOS project is a large redshift survey (Lilly et al. 
2007 ) that was undertaken to study the COSMOS field using 
600 hours of observations with the ESO VLT. COSMOS is an 
HST Treasury Project (Scoville et al. J()07 ) survey of a 2 square 
degree equatorial field with the Advanced Camera for Surveys 
(ACS). It is the largest survey ever completed by HST, utiliz- 
ing 10% (640 orbits) of its observing time over two years (HST 
Cycles 12 and 13), as described further in Koekemoer et al. 
(120071 ). The primary goal of COSMOS is to understand how 
galaxies and AGN evolve over cosmic time in terms of their en- 
vironment, on all scales from groups up to the large scale struc- 
ture of filaments and voids. The COSMOS field is accessible to 
almost all astronomical facilities, which has enabled the com- 
pilation of complete multiwavelength datasets (X-ray, UV, op- 
tical/IR, FIR/submillimeter to radio) to be used in combination 
with the high resolution HST images. 

2.1. Photometric data 

Photometric data in the COSMOS field are available for a wide 
range of wavelengths. In the following, we make use of the 
magnitudes measured in the filters CFHT u* and Ks, Subaru 
Bj, Vj, g^, r^, and z^, and of the Spitzer IRAC magnitudes 
at 3.6 i-im and 4.5 fim. Details about Subaru observations are 
given in Taniguchi et al. ( 120071 ). while data are described 
in McCracken et al. (120091) . Spitzer IRAC data are presented in 
Sanders et al. (^SoDTj). 

Capak et al. (20071 provided a full description of the completed 
data reduction and characteristics of the optical observations. We 
recall that photometry can be optimized by applying offsets to 
the observed magnitudes to reduce the differences between ob- 
served and reference magnitudes computed from a set of tem- 
plate Spectral Energy Distributions (SEDs), as demonstrated by 
Capak et al. (2007; see their Table 13). We adopted the same ap- 
proach, but we derived our own offsets by using the set of SEDs 
that we used to compute absolute magnitudes (see Sect. 13. II ). In 
all cases, the offsets that we derived are similar to those of Capak 
et al. (i2007l ). 

2.2. Morphologies from ACS images 

Morphological parameters for the galaxies are obtained from 
the HST ACS imaging (Koekemoer et al. 12007). The COSMOS 
/ band ACS images have sufficient depth and resolution to al- 
low classical bulge-disk decomposition of L* galaxies at z < 2. 
The size of the COSMOS sample suggests the use of an auto- 
matic and objective morphological classification technique. The 
method adopted for the morphological classification is described 
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in detail in Cassata et al. (|2007l |2009] l and Tasca et al. (12009] ). 
and here we summarize the main steps of the procedure. 
Using a training set of ~ 500 galaxies for which eye-ball mor- 
phological classification is available, the parameters describ- 
ing the galaxy morphology were determined. The classification 
scheme was based on three non-parametric diagnostics of galaxy 
structures (Abraham et al. 2003, Lotz et al. 2004), the concentra- 
tion index C, the asymmetry parameter A and the Gini coefficient 
G, and the galaxy magnitude in the / band (F814W HST/ACS 
band). These parameters were then converted into morphological 
classes: this is done computing in the multi-dimensional parame- 
ter space the distance of each galaxy from objects in the training 
set. The morphological class was then assigned on the basis of 
the most frequent class of the 1 1 nearest neighbours. Given the 
fact that the classification is based on the observed F8 14W band, 
in the redshift range considered in this work (z < 1 ) the effects of 
morphological K-correction are small (Cassata et al. 120071 Tasca 
etal. l2009] l. 

In the following, we use this morphological classification, divid- 
ing galaxies into early-types (including ellipticals and lenticu- 
lars), spirals, and irregulars. 

2.3. Spectroscopic data 

Spectroscopic redshifts in the zCOSMOS survey are obtained 
with the VIMOS spectrograph (Le Fevre et al.'20DT) at the ESO 
Very Large Telescope. The zCOSMOS survey consists of two 
parts. The first part (zCOSMOS-bright) is a pure magnitude- 
limited survey with 15 < / < 22.5 (the / magnitude having been 
measured in the F814W HST/ACS band) and covers the entire 
1 .7 deg^ COSMOS field. This magnitude limit will yield a sam- 
ple of ~ 20 000 galaxies in the redshift range 0.05 i z i 1.2. 
The second part (zCOSMOS-deep) aims to observe ~ 10 000 
galaxies in the redshift range 1.5 i z i 3.0, selected by clearly 
defined color criteria, and is restricted to the central 1 deg^ of the 
COSMOS field. 

Spectroscopic data were reduced with the VIMOS Interactive 
Pipeline Graphical Interface (VIPGI, Scodeggio et al. 120051 ) 
and redshift measurements were derived using the EZ package 
(Garilli et al. 120091) and then visually checked. Each redshift 
measurement was assigned a quality flag, between "0" (impos- 
sible to determine a redshift) and "4" (for which the measure- 
ment is 100% certain); flag "9" indicates spectra with a single 
emission line, for which multiple solutions are possible. Specific 
flags are used to denote Broad Line AGNs. A decimal digit indi- 
cates how closely the redshift agrees with its photometric red- 
shift (Feldmann et al. 2006) computed from optical and near 
IR photometry, using the code ZEBRA (Feldmann et al. 120081 ). 
Further details of the reduction procedure, redshift determina- 
tion, and quality flags are given in Lilly et al. (.2007. .2009,) . 

2.4. Environment 

One of the main scientific objectives of the zCOSMOS survey is 
to study the role of the environment on galaxy evolution to high 
redshift. For the zCOSMOS-bright survey, spectroscopic obser- 
vations were performed with the medium resolution (R ~ 600) 
red grism, which provides a velocity accuracy of ~ 100 km/s. 
Repeated observations of ~ 100 galaxies enabled us to estimate 
the rms velocity uncertainty in each measurement to be of the or- 
der of 1 10 km/s (Lilly et al. 2009 ). The measurement of redshifts 
to such precision allows us to define environments of galaxies 
from the scale of galaxy groups to the larger scales of the cos- 



mic web. 

The density field of the COSMOS survey and the local environ- 
ment of zCOSMOS galaxies have been derived in Kovac et al. 
( |2009al ), where various estimators based on counts in fixed co- 
moving apertures (cylindrical, spherical and Gaussian) and the 
distance to the nearest neighbours are presented, using different 
tracers (flux-limited or volume-limited subsamples). 
In the following, we use overdensities derived by the 5th nearest 
neighbour estimator and computed using volume-limited trac- 
ers. This choice is a good compromise between the smallest ac- 
cessible scales, the robustness of the estimator, and the covered 
redshift range. Bolzonella et al. (2009) discuss in detail the ef- 
fects of different choices, in terms of estimators and tracers, on 
the estimate of the galaxy stellar mass function. 

2.5. Tine zCOSMOS 1 0k briglit sampie 

The zCOSMOS survey is currently ongoing: the data used in the 
present paper are the so-called "10k bright sample" (Lilly et al. 
2009 ), which consists of the first 10 644 observed objects, over 
an area of 1 .402 deg^ with 1 or 2 passes per VIMOS pointing 
and a mean sampling rate of ~ 33%. 

For the present analysis we excluded Broad Line AGNs, stars, 
and objects that had not been included in the statistical sample 
defined in the magnitude range 15 < / < 22.5. We used only 
galaxies with reliable redshifts, i.e. starting from flags 1.5 (see 
Lilly et al. 120071 and ,2009. for details about the flag definition): 
such a sample has an overall reliability of ~ 99%. These objects 
comprise 88% of the overall sample and 95% of the objects in 
the redshift range 0.5 < z < 0.8. 

To obtain a reliable SED fitting when deriving absolute magni- 
tudes and spectrophotometric types (see Sect. 13.1b . we consid- 
ered only objects with apparent magnitudes measured in more 
than 3 photometric bands. 

The final sample used in this paper consists of 8478 galaxies 
satisfying the criteria described above. 



3. Luminosity function estimate 

3. 1 . Absoiute magnitudes and spectropiiotometric types 

Absolute magnitudes were computed following the method de- 
scribed in the Appendix of Ilbert et al. (,2005). The K-correction 
was computed using a set of templates and all available pho- 
tometric information. However, to reduce the template depen- 
dency, the rest-frame absolute magnitude in each band was de- 
rived using the apparent magnitude from the closest observed 
band, shifted to the redshift of the galaxy. With this method, the 
applied K-correction was as small as possible. 
The spectrophotometric types are defined by matching the rest- 
frame magnitudes to the set of templates described in Ilbert et 
al. (2006a}: the four locally observed CWW spectra (Coleman 
et al. 1980) and two starburst SEDs from Kinney et al. ( 119961 ), 
extrapolated toward UV and mid-IR wavelengths, interpolated 
to obtain 62 SEDs, and optimized by using the VVDS spectro- 
scopic data. The SED fitting was performed by minimizing a 
X~ variable on these templates at the spectroscopic redshift of 
each galaxy, providing as output the best-fit spectrum template. 
Galaxies were then classified into four types, corresponding to 
colors of E/SO (type 1), early spirals (type 2), late spirals (type 
3), and irregular and starburst galaxies (type 4). 
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Fig. 1. Evolution of the global luminosity function in the B rest-frame band. Each panel refers to a different redshift bin, which 
is indicated in the label, with also the number of galaxies. The vertical dashed and dotted lines represent the faint absolute limit 
considered in the STY estimate for the zCOSMOS and VVDS sample, respectively. The luminosity functions are estimated with 
different methods (see text for details) but for clarity we plot only the results from (squares) and S TY (solid line). The dotted 
line represents the luminosity function estimated in the redshift range [0.2 - 0.4] and it is reported in each panel as a reference. 
The STY estimates from the VVDS sample are plotted with dashed lines. In the redshift bins [0.8 - 1.0] and [1.0 - 1.3] for the 
zCOSMOS sample we fixed a to the VVDS value. 



3.2. The ALF tool 

Luminosity functions were computed using the "Algorithm for 
Luminosity Function" (ALF), a dedicated tool that implements 
several estimators: the non-parametric l/Vmax (Schmidt 1968), 
(Lynden Bell 119711 in its modified version described by 
Zucca et al. [T997] l. S WML (Efstathiou et al.fTMgjl, and the para- 
metric STY (Sandage, Tammann & Yahil 1979), for which we 
assumed a Schechter function (Schechter il976l . The tool and 
these estimators were described in detail by Ilbert et al. (2005). 
Ilbert et al. (2004) demonstrated that the estimate of the global 
luminosity function can be biased, mainly at the faint end, when 
the measurement band differs in wavelength considerably from 
the rest-frame band in which galaxies are selected. This is caused 
by the fact that, because of the K-corrections, different galaxy 
types are visible in different absolute magnitude ranges at a given 
redshift and fixed apparent magnitude limit. When computing 



the luminosity functions we avoided this bias by using in each 
redshift range (for the C^, SWML, and STY estimates) only 
galaxies within the absolute magnitude range for which the en- 
tire wavelength range of their SEDs is potentially observable. 
We used the complete magnitude range for only the l/Vma.r es- 
timate. Since this estimator underestimates the luminosity func- 
tion for absolute magnitudes fainter than the bias limit (Ilbert et 
al. 120041 1. it provides a lower limit of the faint-end slope. 
Even if this bias is less important when estimating the luminosity 
function of galaxies divided by type, because the K-corrections 
are more similar to each other, we have, however, taken it into 
account. The absolute magnitude limits for the S TY estimate are 
indicated by vertical dashed lines in the figures, and in the tables 
where the best-fit parameters are reported, we provide both the 
total number of objects and the number of galaxies within this 
magnitude limit. 
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Fig. 2. Observed fraction of bright galaxies (Mb < -20.77) of different types as a function of redshift. Error bars are Icr Poisson 
errors. Left panel: spectrophotometric type 1 (red), type 2 (orange), type 3 (green) and type 4 (blue) galaxies. Right panel: morpho- 
logical early-type (ellipticals and lenticulars) (red), spiral (green) and irregular (blue) galaxies; in gray unclassified galaxies. 



3.3. The weighting scheme 

To take account of unknown redshifts (for unobserved objects 
and poor quality spectra), it is necessary to apply a weight to 
each galaxy (Zucca et al. 1 19941 Ilbert et al. 2005). This weight is 
a combination of two different contributions: the target sampling 
rate (TSR) and the spectroscopic success rate (SSR). 
For each galaxy /, we computed the weight w,-, defined to be the 
product of two factors: 

a) wj^^ = 1/TSR = A^phot/A^spec, where A^spec is the number of 
observed sources in the spectroscopic survey and A^phot is the 
number of sources contained in the parent catalogue used to se- 
lect the targets; the TSR thus describes the fraction of sources 
observed in the spectroscopic survey. 

b) wf^ = 1/SSR = N!;IJ(N!;1 - Ntfl), where Nf^l is the 
number of galaxies observed spectroscopically (i.e., excluding 
spectroscopically confirmed stars and Broad Line AGNs) and 
Nlpl^. is the number of objects without a reliable measure of red- 
shift, i.e., "failures". 

Since the galaxies observed spectroscopically were randomly se- 
lected from the parent sample, the TSR is independent of both 
the apparent magnitude and other observed quantities, with an 
approximately constant value of ~ 33%. 

In contrast, the SSR is a function of the selection magnitude, 
which is linked to the signal-to-noise ratio of the spectrum and 
ranges between 97.5% at bright magnitudes and 82% at the 
faintest ones for our subsample of galaxies. Moreover, from the 
photometric redshift distribution is evident that low-confidence 
redshift flags or complete failures correspond to objects at high 
redshift or in redshift ranges where most of the prominent lines 
fall outside the observed wavelength range (Lilly et al. 2009). 
For this reason, we used the Ilbert et al. (2009) release of Zphot 
and computed the SSR in Az - 0.2 redshift bins. Finally, the 
characteristic emission or absorption lines are different for dif- 
ferent galaxy types, as shown in Lilly et al. (2009). We further 
split the computation of SSR in each redshift bin separating red 
and blue galaxies, selected on the basis of their rest-frame U -V 



color We computed the final weights w, — wj x consid- 
ering all the described dependencies. 

This scheme was applied to each observed galaxy, a part from 
two groups of special objects. The spectroscopic catalogue con- 
tains not only randomly targeted objects, because a small frac- 
tion of sources (mainly X-ray sources) were flagged as "com- 
pulsory" when preparing the masks; the TSR of these sources is 
much higher than the global one and was computed separately. 
In some cases, a slit also contained objects in addition to the pri- 
mary target. For these "secondary targets", the SSR was found 
to be far lower than the global relation, because most of these 
sources had not been centered well in the slit, resulting in faint 
spectra. For these objects, the SSR was computed separately. 
The contribution of these two classes of objects to the total sam- 
ple was only ~ 2% and ~ 3%, respectively. 

4. The global luminosity function 

The global luminosity function was computed as a function of 
redshift to z = 1.3, adopting the same redshift bins used in the 
analysis of the VVDS deep field (Ilbert et al. i2005D . to allow a 
direct comparison to be made. In the following, we show the re- 
sults for the Johnson B rest-frame band. 

Figure [T] shows the global luminosity function in redshift bins, 
obtained with the and STY methods. The luminosity func- 
tions derived with the other two methods [l/Vmax and SWML) 
are consistent with those shown in the figures, but are not drawn 
for clarity. The dotted line represents the luminosity function es- 
timated in the redshift range [0.2 - 0.4] and is shown in each 
panel for reference. The S TY estimates from the VVDS sample 
are plotted with dashed lines. Vertical lines represent the bias 
limit described in Sect. 13.21 

The global luminosity function agrees well with the VVDS 
estimates: there are small differences in the normalization in 
some redshift bins, due to the presence of underdense and over- 
dense regions in both samples, which is particularly strong in 
the zCOSMOS sample (see Lilly et al. i2009i ). In particular, the 
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Fig. 3. Luminosity functions of the different spectrophotometric types in redshift bins (indicated in each panel): type 1 in red, type 
2 in orange, type 3+4 in blue, total sample in black. The squares represent the results from the and the solid lines are the results 
from the STY method. The vertical dashed line represents the faint absolute limit considered in the STY estimate. The shaded 
regions represent the 68% uncertainties on the parameters a and M*. 



high zCOSMOS LF normalization in the redshift bin [0.8 - 1.0] 
is due to a prominent structure, discussed in detail in Kovac 
et al. ( |2009ai We note also that the VVDS luminosity func- 
tions were derived with slightly different cosmological param- 
eters (Q.„i - 0.3 and = 0.7). 

Because of the larger area, the zCOSMOS sample is more suit- 
able for constraining the bright part of the luminosity function 
with respect to the VVDS deep sample. With the zCOSMOS 
data, it is also possible to derive some constraints on the M* 
value in the first redshift bin, where for the VVDS sample it was 
necessary to fix M* to the local SDSS value. On the other hand, 
the VVDS fainter magnitude limit allowed us to estimate the 
slope a more reliably at high redshift, where the estimate from 
the zCOSMOS sample is almost unconstrained. For this reason. 



we fixed a to the VVDS value in the redshift bins [0.8 - 1 .0] and 
[1.0-1.3]. 

Having shown full consistency between the zCOSMOS and 
VVDS luminosity functions, we then derived the luminosity 
function in the redshift bins [0.10 - 0.35], [0.35 - 0.55], [0.55 - 
0.75] and [0.75 - 1.00], following the choice of Pozzetti et al. 
(12009 ^) for their study of the zCOSMOS mass function, and in 
Table [T] we report the STY parameters for each bin. 
The errors reported in this table represent a 1 o" confidence level 
for a 1 -parameter estimate {Ax^ - 1.0). These errors underesti- 
mate the true errors. In particular, they are smaller than the pro- 
jection onto the parameter axes of the 2-parameter 68% confi- 
dence ellipses, which, on the other hand, are always an overes- 
timate of the true errors, especially when the errors on the two 
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Fig. 4. Evolution of the luminosity functions for different galaxy types: type 1 (left panel), type 2 (middle panel) and type 3 + 4 
(right panel) galaxies. The colors refer to different redshift ranges: [0.10 - 0.35] in black, [0.35 - 0.55] in cyan, [0.55 - 0.75] in 
magenta and [0.75 - 1.00] in green. The STY estimates are derived with a fixed. The meaning of lines is the same as in Fig. [3] 
points from the estimates are not reported for clarity. 



parameters are highly correlated. 

The luminosity function evolves by ~ 0.7 mag in M* from the 
first ([0.10 - 0.35]) to the last ([0.75 - 1.00]) redshift bins. This 
result agrees with the VVDS results for the same redshift range. 
For the evolution in slope, we find that a is almost constant in 
the first three redshift bins, and then steepens in the last bin. 
To understand how the different galaxy populations contribute to 
this evolution in the global luminosity function, we now quantify 
the contribution of the different galaxy types, by deriving their 
luminosity functions separately. 

5. The contribution of the different galaxy types 

5.1. Spectrophotometric types 

Galaxies were divided into four spectrophotometric classes (see 
Sect. |3.1| above). In the left panel of Fig. |2] we plot the observed 
fraction of bright galaxies for each type as a function of red- 
shift. We identified and selected objects with Mb < -20.77 to be 
galaxies visible for the entire redshift range [0 - 1.0]. The same 
cut was adopted for the VVDS, but in that case it was possible to 
sample the range [0-1.5] because of the fainter magnitude limit 
of the survey. From this figure, it is clear that the bright late-type 
population becomes increasingly dominant at higher redshifts, 
while, correspondingly, the fraction of bright early-type galax- 
ies decreases. These trends in the fraction of the bright galaxy 
population with redshift are similar to those found in the VVDS 
(in the common redshift range), although with some differences 
in the normalizations. The decrease of early-type galaxies in the 
zCOSMOS sample appears less significant than in the VVDS 
sample: starting from similar values at z ~ 0.3, the zCOSMOS 
curve remains significantly (at ~ 2cr level) higher by ~ 20% 
than the VVDS one. This fraction decreases by a factor ~ 2 at 
z ~ 0.9 in zCOSMOS, while it decreases by a factor of ~ 2.3 at 
the same redshift in the VVDS. In the same redshift range, the 
fraction of bright type 3 galaxies increases by a factor ~ 6; the 
strong increase in type 4 galaxies detected in the VVDS is not 
visible here, due to the lower redshift range. 
Luminosity functions were derived for each type in different 
redshift bins. For each type, we derived the luminosity func- 
tion also by fixing a to the value obtained in the redshift range 
[0.30 - 0.80]. This choice allowed us to better constrain the evo- 
lution of M* with redshift and is acceptable because most of 
the a values estimated in the various redshift bins are consistent 



with the [0.30 - 0.80] value. This consistency is marginal for the 
highest redshift bin of type 1 galaxies, but in this bin the STY 
estimate is poorly constrained because the faint end of the lumi- 
nosity function is inadequately sampled, due to the magnitude 
limit of the survey. 

For type 1, 2 and 3 galaxies, this value of a agrees with the 
VVDS estimate. For type 4 galaxies, the parameter a is uncon- 
strained for z > 0.3: for this reason, we fixed a to be the value 
derived from the first redshift bin, which is consistent with the 
VVDS value for type 4 galaxies. 

The evolution of M* and 0* with redshift for the different types 
are consistent with the VVDS results in the common redshift 
range. 

Given the uncertainties in the luminosity function estimate of 
type 4 galaxies, we repeated the analysis by grouping together 
type 3 and 4 galaxies: the results are shown in Fig. [3] and the 
S TY parameters for each bin are reported in Table[T] for both a 
free and a fixed. Parameters for the reference bin [0.30 - 0.80] 
are also reported in the table. In Fig.[3]the squares represent the 
results from the and the solid lines are the results from the 
STY method: type 1 galaxies are shown in red, type 2 galaxies 
in orange, type 3+4 galaxies in blue. The total sample is shown 
in black. The shaded regions represent the 68% uncertainties in 
the parameters a and M*. From this figure, it is clear that at low 
redshifts (z < 0.35) late-type galaxies dominate for faint mag- 
nitudes {Mb > -20), while the bright end is populated mainly 
by type 1 galaxies. At higher redshift, late-type galaxies evolve 
strongly and at redshift z > 0.55 the contributions of the various 
types to the bright end of the luminosity function are compara- 
ble. The faint end remains dominated by late-type galaxies over 
the entire redshift range. 

To visualize the evolution with redshift, in Fig. |4] we plot the 
luminosity functions in different redshift bins for each galaxy 
type: type 1 galaxies in the left panel, type 2 galaxies in the mid- 
dle panel, and type 3+4 galaxies in the right panel. The different 
colors represent different redshift bins: [0.10 - 0.35] in black, 
[0.35-0.55] in cyan, [0.55-0.75] in magenta, and [0.75 - 1.00] 
in green. To follow the evolution in M* and (p*, we show the 
STY estimates obtained with a fixed. For type 1 galaxies, evo- 
lution occurs in both luminosity and normalization: M* bright- 
ens by ~ 0.6 mag and (/>* decreases by a factor ~ 1.7 between 
the first and the last redshift bin. Type 3+4 galaxies also evolve 
both in luminosity and normalization, but with an opposite trend 
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Fig. 5. Luminosity functions of the different morphological types in redshift bins: early types (ellipticals and lenticulars) in red, 
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for the normalization: a brightening with redshift of ~ 0.5 mag 
is evident in M*, while <p* increases by a factor ~ 1.8. Type 2 
galaxies exhibit a milder evolution, involving a brightening of 
~ 0.25 mag in M* and no significant evolution in (p*. 
The galaxy stellar mass functions of the various types (Pozzetti 
et al. |20091 l show differences in the massive part stronger that the 
differences we find in the bright part of the luminosity functions. 
This is due to the fact that the mass-to-light ratio of early-type 
galaxies is on average higher than that of late-type galaxies: in 
fact, the contribution to the bright end of our B band luminosity 
functions comes from both massive red galaxies and blue galax- 
ies with strong star formation. This fact also implies that the bi- 
modality observed in the global galaxy stellar mass function by 
Pozzetti et al. ( I2009I I is not detected in the global luminosity 
function. 



5.2. Morphological types 

A major advantage of the zCOSMOS survey is the availabil- 
ity of galaxy morphologies obtained from the HST ACS im- 
ages (Koekemoer et al. 2007). Galaxies were divided into early 
types (including ellipticals and lenticulars), spirals, and irregu- 
lars following the classification described in Sect. l2.2l above. In 
the right panel of Fig. |2] we plot the observed fraction of bright 
galaxies of each morphological type as a function of redshift. 
Unclassified galaxies are a very small fraction of the sample 
(~ 3%) and are uniformly distributed with redshift. 
As observed for spectrophotometric types, late-type galaxies 
steadily increase their fraction with increasing redshift, while 
the fraction of early types decreases. However, at low redshift 
the fraction of morphologically classified early types is higher 
than that of type 1 galaxies, as is evident by comparing the left 
and right panels in Fig. |2] 
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The luminosity functions of the different types are shown in Fig. 
|5] (early types in red, spirals in green, irregulars in blue) and the 
STY parameters for each bin are reported in Table [T] where we 
show results for both a free and a fixed to the value determined 
in the redshift range [0.30 - 0.80]. 

At low redshift (z < 0.35), early-type galaxies dominate the 
bright end of the luminosity function, while spiral galaxies domi- 
nate the faint end. Irregular galaxies increase their contribution at 
the lowest luminosities. At intermediate redshift ([0.35 - 0.75]), 
spiral galaxies increase their luminosities and their contribution 
to the bright end of the luminosity function is similar to that of 
the early types. At high redshift (z > 0.75), irregular galaxies 
evolve strongly and the three morphological types contribute al- 
most equally to the total luminosity function. Irregular galaxies 
show an evolution of a factor ~ 3.3 in 0* from low to high red- 
shift. This evolution occurs mainly in the last redshift bin, while 
for z < 0.75, the contribution of these galaxies to the global lu- 
minosity function is significantly lower than that of spirals and 
early types. 

Scarlata et al. (2007') derived luminosity functions for COSMOS 
galaxies, using photometric redshifts and the ZEST morpholog- 
ical classification. The shape parameters we find for early types 
and spirals are consistent with those found by Scarlata et al. 
(I2007I I for their early-type and disk galaxies, but are significantly 
different for irregular galaxies, which have a much flatter slope 
in Scarlata et al. (2007). These diff'erences are likely due to the 
different morphological classification applied. 
Although the general trend in the luminosity functions of the 
different morphological types is similar to the results obtained 
in the previous section for spectrophotometric types, some dif- 
ferences are present. In particular, there are more morphological 
early-type than type 1 galaxies at the faint end of the luminos- 
ity function (Mb xL - 19.5 in the first redshift bin). In the fol- 
lowing, we discuss the relationship between spectrophotometric 
and morphological types, paying particular attention to early- 
type galaxies. 

6. Spectrophotometric versus morphological types 

Although there is a broad agreement between our results for 
spectrophotometric and morphological types, it is reasonable 
that there should also be some differences. Spectrophotometric 
types are based on the galaxy SEDs and therefore depend on the 
star formation history, while morphological types reflect mainly 
the dynamical history of the galaxy. These classifications can 
also be affected by different observational and methodological 
biases. 

We compared the spectrophotometric and morphological results 
for each galaxy. Considering early-type galaxies, we have 2387 
morphological early types, 1504 of which (63%) are classified 
as type 1 galaxies. In contrast, 71% of the 2105 type 1 galaxies 
are classified as morphological early types. If we consider only 
galaxies whose SED is most accurately reproduced by the most 
extreme (i.e. reddest) type 1 template, the fraction increases to 
80%. The remaining fraction of type 1 galaxies that are not clas- 
sified as morphological early types is in part due to objects being 
without morphological classification (~ 3%) and in part due to 
a population of "red" spirals, many of which on visual inspec- 
tion appear to be edge-on spiral galaxies, often dominated by a 
strong dust lane (see Tasca et al. 2009 for a detailed discussion). 
The red SED for these galaxies is probably caused by a signifi- 
cant amount of dust extinction. 

We then considered in more detail the ~ 37% of morphological 
early-type galaxies that were not classified as spectrophotomet- 



ric type 1 . We termed "blue" early types those with spectropho- 
tometric type 2, 3 or 4; among the morphological early types 
with spectrophotometric type 1, we termed "very red" those best 
fitted by the reddest type 1 template and "red" those best fit- 
ted by the other type 1 templates. These "blue" early types were 
visually inspected (see details in Tasca et al. |2009l l revealing a 
class of face-on late-type galaxies with morphological parame- 
ters typical of an early-type population. 

To explore the properties of these classes of morphological early 
types, composite spectra were generated for these objects (fol- 
lowing Mignoli et al. I2009I I by averaging their spectra and di- 
viding them into bins brighter and fainter than Mb = -20.77. 
These spectra are shown in Fig. |6] where the Y-axis is rescaled 
arbitrary for clarity, for the bright sample on the left and for 
the faint sample on the right. Blue, magenta, and red spectra 
correspond, respectively, to "blue", "red", and "very red" early 
types. From this figure, clear differences are visible in the spec- 
tra, in particular prominent emission lines are present in "blue" 
early types. By comparing the bright and the faint "blue" early 
types, we find that the Ha line equivalent width increases for 
fainter galaxies and that the ratio between [Nil] and Ha lines is 
stronger for brighter galaxies. The [NH] / Ha and [OHI] / Hp 
ratios of the bright "blue" early types are consistent with those 
of liners. Bright "blue" early types are ~ 30% of the total num- 
ber of bright early types; by considering the faint early types, the 
fraction of "blue" early types increases to ~ 44%. 
Similar spectra for "blue" early types at low (log M/Mq < 9) 
and high (log Ml Mq > 10.5) masses are shown by Pozzetti et 
al. (2009), who find a more evident separation between the prop- 
erties of objects in these two subclasses. 

As a final check we also considered a sample of bona-fide early- 
type (ET) galaxies (Moresco et al. I2009I I selected in a more 
conservative way by combining information on morphologies, 
spectrophotometric types, colors, and emission line equivalent 
widths (see Pozzetti et al. 2009 for the exact criteria adopted). 
This conservative selection reduces the number of objects in the 
early-type sample to 98 1 : therefore, we use bona-fide ET galax- 
ies mainly as a comparison sample. 

We now explore the luminosity function of early-type galaxies. 

6. 1 . The luminosity function of eariy-type galaxies 

In Fig.|2]we show the luminosity function of early-type galaxies 
from different samples: bona-fide ET (in black), type 1 (in red), 
and morphological early-type (in green) galaxies. 
The most apparent feature in this figure is the excess of faint 
(Mb i - 19.5) morphological early-type galaxies at z < 0.55. 
This effect is not evident at higher redshift because this faint pop- 
ulation is cut by the magnitude limit of the survey. This excess is 
due to the population of "blue" early types discussed in the pre- 
vious section. A similar effect was already noticed by Ilbert et al. 
(2006b), who found that the main contribution to the faint end 
of the luminosity function of bulge-dominated galaxies comes 
from blue objects. 

Moreover, it is interesting that the bright ends of the luminosity 
functions of type 1 and morphological early-type galaxies are al- 
most indistinguishable at all redshifts. 

We also note that type 1 and bona-fide ET galaxies luminos- 
ity functions have similar shapes, but different normalizations. 
The slopes are consistent to within the errors, as well as the M* 
values; the difference in (p* is about a factor of 2, due to the se- 
lection criteria adopted in constructing the bona-fide ET sample. 
The similarity between these luminosity functions implies that, 
even if the type 1 galaxy sample is contaminated by a fraction 
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Fig. 8. Left column: luminosity functions in the lowest (blue) and highest (red) quartile of the density distribution in redshift bins 
for the total sample. Overdensities are derived from the 5 nearest neighbours. Middle and right columns: luminosity functions in 
the lowest (middle column) and highest (right column) quartile of the density distribution in redshift bins. In each panel the global 
luminosity function is shown in black, while the luminosity function of type 1 and type 3+4 galaxies are drawn in red and blue, 
respectively. The meaning of symbols and lines is the same as in Fig. [3] The shaded regions represent the 68% uncertainties on the 
parameters a and M*. For type 1 and type 3+4 samples the S TY estimates with a fixed are shown and therefore the shaded area is 
not drawn. 



of "red" spirals, the contribution of this population is approxi- 
mately constant with both luminosity and redshift. 

7. The role of the environment 

To understand the effect of the environment, we derived the lu- 
minosity functions for galaxies in underdense and overdense re- 
gions, using the density estimates described in Sect. 12.21 We re- 
peated the analysis using various density estimators, finding that 
the general trends are similar in all cases. In the following, we 
show the results obtained for overdensities derived with the 5th 
nearest neighbour estimator 



To emphasize the differences between the underdense and over- 
dense environments, we consider only the two extreme tails of 
the overdensity distribution. We also ensure that we had roughly 
the same number of objects in each subsample, to avoid spuri- 
ous normalization effects. For this reason, we divided the galaxy 
density distribution into quartiles, within each redshift bin, and 
derived the luminosity functions for each subsample. In the fol- 
lowing, we compare the results for the lowest and highest quar- 
tiles of the overdensity distributions. The use of both quartiles 
and subsamples reduces the number of usable objects, and we 
therefore decided to use only three redshift bins, [0.1 - 0.4], 
[0.4 - 0.7], and [0.7- 1.0]. 
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Fig. 9. Evolution of the luminosity functions for different galaxy types in different environments (low density on the left, high 
density on the right): type 1 (upper panels) and type 3+4 (lower panels) galaxy samples. The color refer to different redshift ranges: 
[0.10 - 0.40] in black, [0.40 - 0.70] in blue and [0.70 - 1.00] in green. The S TY estimates are derived with a fixed. The meaning 
of lines is the same as in Fig. [3] points from the estimates are not reported for clarity. 



The most appropriate choice of quartile values depends on the 
topic being investigated. To measure the contribution of differ- 
ent galaxy types to the global luminosity function in different 
environments, it is necessary to derive the quartile values for 
the total sample and use the same density cuts for all the sub- 
samples. However, if the aim is to compare the behaviour of the 
luminosity function of specific subsamples of galaxies in under- 
dense and overdense environments, we must carefully select for 
each subsample the appropriate quartiles (i.e., as derived for that 
particular subsample). 

The use of quartiles makes difficult to derive a number den- 
sity from the (p* parameter of the luminosity function, because 
this kind of selection does not conserve information about the 
volume occupied by the overdensity associated to each galaxy. 
Therefore, when analysing luminosity functions in different en- 
vironments, only a shape comparison is allowed. However, the 
comparison of 0* values is correct when luminosity functions for 
galaxies in the same quartile are compared. 
In the following, we show results both for the total sample and 
for early- and late-type galaxies, using the type 1 and type 3-h4 
subsamples. For the early-type galaxies we decided to use the 
type 1 instead of the bona-fide ET sample for three reasons: 1 . 
the more reliable statistical analysis possible for a higher num- 
ber of objects; 2. the homogeneity in the classification compared 
to the late-type galaxy sample; 3. the similarity between the type 
1 and bona-fide ET galaxy luminosity functions (see previous 
section). However, we checked that the results obtained with 
the bona-fide ET sample (although with larger uncertainties) are 
consistent with those derived in the following for the type 1 sam- 



ple. 

In the first column of Fig. [8] we show the luminosity function for 
the total sample in the lowest (blue lines and points) and high- 
est (red lines and points) quartiles of the density distribution in 
redshift bins. The luminosity function of galaxies in overdense 
regions is consistently of brighter M* : the slope in underdense 
regions is steeper in the first and second redshift bins, while the 
a values in the highest redshift bin, in both environments, are 
consistent with each other, in the significant margins of large er- 
ror Increasing the redshift, the contribution from the luminosity 
function of galaxies in underdense environments increases and 
the crossing point between the two LFs progressively moves to- 
wards brighter magnitudes. To explore whether this behaviour 
is due to the different morphological mixes in the two environ- 
ments, in the second and third column of Fig. [8] we show the 
luminosity function in underdense and overdense environments: 
the total luminosity function is shown in black, while the contri- 
bution from type 1 and type 3+4 galaxies are drawn in red and 
blue, respectively. In this case, the quartile values are those de- 
rived for the total sample. The parameters of the S TY estimates 
are reported in Table |2] The STY parameters for type 1 and type 
3+4 galaxies samples, although formally constrained, are poorly 
determined. We therefore repeated the STY estimate by fixing 
a for each type to the global value derived in Sect. 15.11 after 
verifying that this choice is acceptable. In low density environ- 
ments, the main contribution to the luminosity function is from 
type 3+4 galaxies, while for high density environments an im- 
portant contribution is that of type 1 galaxies at the bright end. 
The differences between the global luminosity function in the 
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two environments are due not only to the different relative num- 
bers of type 1 and type 3+4 galaxies, but also to their relative 
luminosity distributions. The value of M* in underdense regions 
is always fainter than in overdense environments, by ~ 0.50 mag 
and ~ 0.65 mag for type 1 and type 3+4 galaxies, respectively. 
For 0*, between underdense and overdense environments there 
is a decrease of a factor ~ 2 for type 3+4 galaxies and an in- 
crease of ~ 2 for type 1 galaxies (except in the last redshift bin, 
where the S TY estimate is poorly determined). These results in- 
dicate that galaxies of the same type in different environments 
have different properties. 

We now investigate if there is also a differential evolution for dif- 
ferent environments within each class of objects. To achieve this, 
we need to use, for each subsample, the appropriate quartiles 
(i.e., those derived for that particular subsample). The parame- 
ters of the STY estimates obtained with this choice are reported 
in Tables 

To visualize the evolutionary effects more clearly, in Fig. |9] we 
plot in the same panel the luminosity functions in different red- 
shift bins for each galaxy type in the two environments, showing 
for clarity only the STY estimates. For type 3+4 galaxies, the 
evolution in the luminosity function within underdense and over- 
dense regions is similar for the first and the second redshift bin, 
and occurs mainly in luminosity: there is a brightening of ~ 0.5 
mag in M* for all cases, and a slight variation in 0*. In contrast, 
between the second and the third redshift bins there is little lumi- 
nosity evolution, but a significant evolution in <p*, which differs 
in the two environments: in overdense regions 0* increases by a 
factor ~ 1 .6, while in underdense regions this increase reaches a 
factor -2.8. 

By analyzing the type 3+4 galaxy population in more detail, 
we find that the evolution in the relative ratio of type 3 to type 
4 galaxies differs between the two environments: in particu- 
lar, moving from the second to the third redshift bin, the ratio 
Ntype4/N,ype^ for bright {Mb < -20.25) galaxies changes from 
33% to 46% in low density environments, but remains almost 
constant (~ 35%) in high density environments. Therefore, the 
strong evolution observed for type 3+4 galaxies from z ~ 0.5 to 
z ~ 0.9 is mainly due to type 4 galaxies. The adopted limit in ab- 
solute magnitude corresponds to the bias limit (see Sect. 13. 2t in 
the highest redshift bin and allows a direct comparison of galaxy 
numbers at different redshifts. 

For type 1 galaxies, the situation is more unclear, because the 
lower number of objects introduces larger uncertainties in the 
luminosity function estimates. Between the first and the second 
redshift bin, there is a similar evolution in both environments, 
while passing from the second to the third redshift bins there 
are indications of a more significant evolution in underdense re- 
gions. 

To check if the choice of different quartiles for each redshift bin 
can affect these results, we repeated the analysis using in each 
redshift bin the value of the quartiles derived for the whole red- 
shift range [0.1 - 1], finding that the trends do not change. 
These results agree with those found studying the galaxy stel- 
lar mass function in different environments (Bolzonella et al. 
12009 ). In particular, in high density environments there is a 
strong contribution of early-type galaxies, which are dominant 
at high masses in the galaxy stellar mass function and at bright 
magnitudes in the luminosity function: however, this dominance 
is clear at all redshifts in the galaxy stellar mass function, 
but decreases with increasing redshifts in the luminosity func- 
tion. This is due to the fact that at high redshift blue galax- 
ies brighten because they increase their star formation activity. 
Finally, Bolzonella et al. (.2009J find that the contribution of 



early-type galaxies to the total galaxy stellar mass function in- 
creases more rapidly in high density environments: consistently, 
we find that the relative contribution of type 1 and type 3+4 
galaxies to the bright end of the luminosity function is similar 
at high redshift in both environments, but at low redshift type 1 
galaxies dominate mainly at high density. 

8. Summary and conclusions 

An unbiased and detailed characterization of the luminosity 
function is the basic requirement in many astrophysical issues 
and it is of particular interest to assess the role of environment 
on the evolution in the luminosity function of different galaxy 
types. With this aim, we have studied the evolution in the lumi- 
nosity function to redshift z ~ 1 for the zCOSMOS 10k sample, 
for which both accurate galaxy classifications and a detailed de- 
scription of the density field are available. 

The main results of this analysis are the following: 

- Global luminosity function: The global luminosity function 
shows a brightening of ~ 0.7 mag in M* from z ~ 0.2 to z ~ 0.9, 
in agreement with the VVDS results in similar redshift ranges. 

- Luminosity functions by spectrophotometric types: To quantify 
the contribution of the different spectrophotometric types, we 
have considered three classes of objects: type 1 galaxies, corre- 
sponding to early-type SEDs, type 3+4 galaxies, corresponding 
to late-type SEDs, and the intermediate class of type 2 galax- 
ies. At low redshift (z < 0.35), type 3+4 galaxies dominate the 
luminosity function at faint magnitudes {Mb > -20), while the 
bright end is populated mainly by type 1 galaxies. At higher red- 
shift, type 3+4 galaxies evolve strongly and therefore, at redshift 
z ~ 1, the contribution to the bright end of the luminosity func- 
tion of the various types is comparable. The faint end remains 
dominated by type 3+4 galaxies over the entire redshift range. 
For type 1 galaxies, the evolution occurs in both luminosity and 
normalization: not only does M* brighten by ~ 0.6 mag but <p* 
also decreases by a factor ~ 1.7 between the first and the last 
redshift bin. Type 3+4 galaxies also evolve in both luminosity 
and normalization, but with an opposite trend for the normaliza- 
tion: a brightening with redshift by ~ 0.5 mag is present in M*, 
while 0* increases by a factor ~ 1.8. Type 2 galaxies show a 
milder evolution, of a brightening of ~ 0.25 mag in M* and no 
significant evolution in (p*. 

- Luminosity functions by morphological types: At low red- 
shift (z < 0.35), morphological early-type galaxies dominate the 
bright end of the luminosity function, while spiral galaxies domi- 
nate the faint end. Irregular galaxies increase their contribution at 
the lowest luminosities. At intermediate redshift ([0.35 - 0.75]), 
spiral galaxies increase their luminosities and their contribution 
to the bright end of the luminosity function is similar to that of 
the morphological early-type one. At high redshift (z > 0.75), ir- 
regular galaxies evolve strongly and therefore the three morpho- 
logical types contribute almost equally to the total luminosity 
function. Irregular galaxies exhibit an evolution of a factor ~ 3.3 
in (p* from low to high redshift: this evolution occurs mainly in 
the last redshift bin, while, for z < 0.75, the contribution of these 
galaxies to the global luminosity function is significantly lower 
than that of spirals and early types. 

- The role of the environment: For the total sample, the lumi- 
nosity function of galaxies in overdense regions always has a 
brighter M* and a flatter slope. Increasing the redshift, the con- 
tribution to the luminosity function of galaxies in underdense 
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environments increases, and the crossing point between the two 
LFs moves progressively towards brighter magnitudes. In low 
density environments, the main contribution to the luminosity 
function originates in type 3+4 galaxies, while for high den- 
sity environments there is an important contribution from type 
1 galaxies to the bright end. 

The differences between the global luminosity functions in the 
two environments are due not only to a difference in relative 
numbers of type 1 and type 3+4 galaxies, but also to their relative 
luminosity distributions. The value of M* in underdense regions 
is always fainter than in overdense environments, by ~ 0.50 mag 
and ~ 0.65 mag for type 1 and type 3+4 galaxies, respectively. 
For (p*, between underdense and overdense environments, there 
is a decrease of a factor ~ 2 for type 3+4 galaxies and an increase 
of ~ 2 for type 1 galaxies. These results indicate that galaxies of 
the same type in different environments have different proper- 
ties. 

We have also detected differential evolution for type 3+4 galax- 
ies in different environments. We found that the evolution in their 
luminosity function in underdense and overdense regions is sim- 
ilar between z ~ 0.25 and z ~ 0.55, and is mainly in luminosity: 
there is a brightening of ~ 0.5 mag in M* in all cases, but only a 
slight variation in (p*. In contrast, between z ~ 0.55 and z ~ 0.85 
there is little luminosity evolution but a strong evolution in <p*, 
which differs in the two environments: in overdense regions, 0* 
increases by a factor ~ 1 .6, while in underdense regions there is 
an increase of a factor ~ 2.8. 

Analyzing the type 3+4 galaxy population in more detail, we 
find that the evolution in the relative ratio of type 3 to type 4 
galaxies differs between the two environments: in particular, be- 
tween z ~ 0.55 and z ~ 0.85, the ratio Ntype4/Niypei for bright 
(Mb < -20.25) galaxies changes from 33% to 46% in low den- 
sity environments, but remains almost constant (~ 35%) in high 
density environments. Therefore, the strong evolution observed 
for type 3+4 galaxies between z ~ 0.55 and z ~ 0.85 is mainly 
due to type 4 galaxies. 

For type 1 galaxies, the situation is unclear, because the lower 
number of objects introduce larger uncertainties in the luminos- 
ity function estimates. Between the first and the second redshift 
bins, there is a similar evolution in both environments, while 
passing from the second to the third redshift bins, there are indi- 
cations of a more significant evolution in underdense regions. 



All these results are summarized in Fig.[TOl where the evolution 
of the parameters M* and (p*, obtained with a fixed, is shown for 
different samples; the label near each point indicates the mean 
redshift of the considered bin. The upper panel of Fig.fTOlrefers 
to galaxies of different spectrophotometric types in the total sam- 
ple: type 1 galaxies in red, type 2 galaxies in orange and type 3+4 
galaxies in blue. From this figure it is immediately evident the 
symmetric behaviour of type 1 and type 3+4 galaxies, as well as 
the almost negligible evolution of type 2 galaxies. While we find 
that the M* evolution for type 1 and type 3+4 galaxies is con- 
sistent with passive evolution, it is more interesting the density 
evolution. Given the fact that (p* is proportional to the density of 
M* galaxies, this figure is consistent with a scenario where a part 
of type 3+4 galaxies is transformed in type 1 galaxies with in- 
creasing cosmic time, without significant changes in the fraction 
of intermediate-type galaxies. 

The middle and lower panels of Fig.[TO]show the evolution of the 
parameters M* and (p* for galaxies of different types in different 
environments: type 1 galaxies in underdense (magenta triangles) 
and overdense (red circles) environments in the middle panel. 
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Fig. 10. Evolution of the M* - (p* parameters (derived with a 
fixed) for different samples: the label near each point indicates 
the mean redshift of the considered bin. Upper panel: type 1 
(red), type 2 (orange) and type 3+4 (blue) galaxies in the total 
sample. Middle panel: type 1 galaxies in overdense (red circles) 
and underdense (magenta triangles) environments. Lower panel: 
type 3+4 galaxies in overdense (blue circles) and underdense 
(cyan triangles) environments. 
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type 3+4 galaxies in underdense (cyan triangles) and overdense 
(blue circles) environments in the lower panel. The behaviour of 
type 1 and type 3+4 galaxies in overdense environments is con- 
sistent with the transformation scenario described above for the 
total sample. In underdense environments, a similar scenario is 
acceptable, with the exception of the highest redshift bin, where 
type 1 galaxies show a strong increase of ^* with respect to the 
previous bin. 

The lower panel of Fig.[TO]clearly shows the different evolution 
which type 3+4 galaxies are undergoing in the different envi- 
ronments, with a much stronger density evolution in underdense 
regions. 

This indicates that the bulk of the tranformation described above 
from blue to red galaxies in overdense regions probably hap- 
pened before z ~ 1 , while it is still ongoing at lower redshifts in 
underdense environments. 
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Table 1. S TY parameters (with Icr errors) in different redshift bins for total sample and different galaxy type subsamples. 
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Table 2. STY parameters (with Icr errors) for different galaxy samples in different envirormients: total sample quartiles. 

sample z-bin Number<°) Number"'' a M^^^ - 51og(/i7o) ^'ClO'^ /i^p Mpc"^) 
Under-dense environments 

Total sample 0.10-0.40 680 655 ^TOAA^^ 1.62+;j-^^ 

Type 1 galaxies 0.10-0.40 92 92 -0.2S.+^H -20.34+2-39 0.59*ff^ 

-0.30 fixed -20.37+0-18 ^-^^^w. 

Type 3 -I- 4 galaxies 0.10-0.40 499 477 -1.33+o [2 -20.02+o-| ^-^^-bu 

-1.47 fixed -20.3r|}| 0.85+g'g^ 

Total sample 0.40-0.70 675 628 ^TW^^ -20.97+^^1 0.96+" " 

Type 1 galaxies 0.40-0.70 131 127 -0.88!°-^^ -21.19!?-| 0.21-^,^^ 



0.20 ^"-^'-O.SS "-^"-0.32 

+0.47 Ti 1Q+0.62 n 91 +0.12 

-0.30 fixed -20-50!21^ ^-^"^^oqI 

Type 3 -I- 4 galaxies 0.40-0.70 449 413 -1.64+°|« -21-26+™ 0-3''l l 

-1.47 fixed -20.98+°-}^ 0.55+i;-°^ 

Total sample 0.70-1.00 606 594 -0.93+"-^^ -20.73+^ | 2.82+|| 

Type 1 galaxies 0.70-1.00 129 129 0.38+of -20.34+02' 0.54+|||^ 

-0.30 fixed -20.64+o:|i 0.64+g«« 

Type 3 -I- 4 galaxies 0.70-1.00 365 357 -1.51+2*^ "20-93^^2 ^■^'^m 

-1.47 fixed -20.90;ii-|° 1.40;°f 

Over-dense environments 

Total sample 0.10-0.40 713 678 -0.86!" "^ -20.90!" " 2.14+"-^' 

Type 1 galaxies 0.10-0.40 264 262 -0.12+°-|| -20-50+I ^ 1-63 

-0.30 fixed -20.72:|^ 1-44 fioo 

Type 3 + 4 galaxies 0.10-0.40 327 298 -1-34!2 |^ -20.78+°-^ O.Sl^JJ I^ 

-1.47 fixed -21.10!°-| 0.33+°-°^ 

Total sample 0.40-0.70 703 684 -0.81!" (^ -21.14!|j-|^ 1.34+|| 

Type 1 galaxies 0.40-0.70 263 256 -0.59+||f -21-36+[li 0.49+^5^ 

-0.30 fixed -21.05+™ O-^Oqm 

Type 3 -I- 4 galaxies 0.40-0.70 306 298 -1 •04+0-29 -20-93+!!-?^ O-^'^p i 

-1.47 fixed -21.49;°-|^ 0.32+|°| 

Total sample 0.70-1.00 656 636 -1.01!"-| -21.44!|-|^ 1.44+^-| 

Type 1 galaxies 0.70-1.00 207 202 -0.48+ff» "21-35!gi 0-49!of 

-0.30 fixed -21.24+™° 0.51+||-|j* 

Type 3 -I- 4 galaxies 0.70-1.00 300 288 -1.34+«-« -21-38+!!-2 ^-^Kll 

-1.47 fixed -21.49!°]^ 0.58!°-| 

(a) Number of galaxies in the redshift bin 

(b) Number of galaxies brighter than the bias limit (sample used for 5 TY estimate; see the text for details) 
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Table 3. S TY parameters (with Icr errors) as a function of the environment (from the 5 nearest neighbours) in different redshift 
bins, for different galaxy types. In each sample the appropriate quartiles are used. 



z-bin quartile Number^"' Number**^ 



or M;^„ - 51og(/i7o) K Mpc-^) 



Total sample 



-1 16+" 

-!!QS 
-0.86+!o' 



-20 44+"-" 

-m 



T62^ 
2.14 



0.10-0.40 first 680 655 

0.10-0.40 fourth 713 678 



-20.90 



m- 



0.40-0.70 first 675 628 -1.20! 

0.40-0.70 fourth 703 684 -0.81 



-20.97: 
-21.14+P^ 



0.96 
1.34 



0.70-1.00 first 606 594 -0.93 

0.70-1.00 fourth 656 636 -1.01 



-20.73 
-21.44^ 



■n.is 



2.82_ 

1 44+0.24 



Type 1 galaxies 



+0.20 
-0.19 

-0.30 fixed 
-0 20+°"^ 
-0.30 fixed 



-20 34+"-^'' 
-20 25+" " 
-20 65+°-^' 
-20 79+"- '2 

'^-0.12 



+0.14 

-A& 

-0.08 
Q gj+aX)9 

90+^-''''^ 



0.10-0.40 first 160 160 

0.10-0.40 fourth 166 164 



-0.37 



.0 39+"-^<» 

-0.30 fixed 
.0 99+0.28 

"•^^-0.27 

-0.30 fixed 



-20.65-- 
-20.58!»:«' 
-21 98+''-50 

-21.05+«;? 



0.96 
1.01 



"•^-'-0.08 

45+"-''^ 
20+'' '° 

-QiQ 
42+°*^ 
74+0- 14 

"•'^-0.20 

0.76+!;!!« 



0.40-0.70 first 183 179 

0.40-0.70 fourth 187 180 



-0 04*" " 
-0.30 fixed 

-0.30 fixed 



0.70-1.00 first 160 160 

0.70-1.00 fourth 166 163 



-20 55+0-30 
-20.68!0[? 

-21 23+™ 



0.39 
0.42 



7m 



Type 3+4 galaxies 



— 1 O9+0.U 
-1.JZ._QJ3 

-1.47 fixed 
-1 34+0" 

-0.11 

-1.47 fixed 



a+0.23 

-20.28:11 
-20.64:|| 
-20 94+017 

^"•^ -0.19 

-21 40+"'* 
-20.96+01^ 
-20.92!»:« 



'^•'^^-0.32 

14^°°'^ 
0.78+1! 

0-52tii] 
-023+03S- 

0:46+|{^ 

0.84!»| 
41+0.02 

"•^^-0.02 



0.10-0.40 first 422 406 

0.10-0.40 fourth 439 407 



-19.98+ 



— 1 7^+0.28 

-0 ^7 

-1.47 fixed 
-1 04+0 25 

-1.47 fixed 



0.40-0.70 first 375 343 

0.40-0.70 fourth 391 382 



5+0.52 
-0.50 

-1.47 fixed 

.1 44+0.38 
■ -0.38 

-1.47 fixed 



-21 OO+o^ 

-20.89+o|o 
-21.49+0:^2 



0.70-1.00 first 328 324 

0.70-1.00 fourth 357 342 



-1.63+ 



-21.52+: 



1.30 
0.67 
0.65+: 
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(a) Number of galaxies in the redshift bin 

(b) Number of galaxies brighter than the bias limit (sample used for 5 TY estimate; see the text for details) 



